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The quantum statistics of damped optical solitons is studied using cumulant-expansion tech- 
niques. The effect of absorption is described in terms of ordinary Markovian relaxation theory, by 
coupling the optical field to a continuum of reservoir modes. After introduction of local bosonic 
field operators and spatial discretization pseudo-Fokker-Planck equations for multidimensional s- 
parametrized phase-space functions are derived. These partial differential equations are equivalent 
to an infinite set of ordinary differential equations for the cumulants of the phase-space functions. 
Introducing an appropriate truncation condition, the resulting finite set of cumulant evolution equa- 
tions can be solved numerically. Solutions are presented in a Gaussian approximation and the quan- 
tum noise is calculated, with special emphasis on squeezing and the recently measured spectral 
photon-number correlations [Spatter et at, Phys. Rev. Lett. 81, 786 (1998)]. 



I. INTRODUCTION 



The study of the nonclassical properties of optical 
pulses propagating in linear and nonlinear media has 
been a subject of increasing interest In particu- 

lar effects such as quadrature-phase squeezing ]l],f§-f), 
photon- number squeezing [f^-|To|], intraspectral quantum 
correlations jy], and entanglement JlJ,[l3| offer novel 
possibilities for optical communication and quantum 
nondemolition (QND) measurement [§],[§. In this connec- 
tion the development of reliable sources of nonclassical 
light has been of particular interest. Promising candi- 
dates for such sources have been nonlinear optical ele- 
ments converting coherent light produced by lasers into 
nonclassical light. 

Since nonclassical properties are always in competi- 
tion with the noise associated with unavoidable losses, a 
quantum theoretical description is required that neces- 
sarily includes absorption. Numerous approaches to the 
problem of quantization of radiation in absorbing media 
have been made (see, e.g., |Q and references therein). 
In most quantization schemes the influence of absorp- 
tion has been omitted. Concepts have been developed for 
both dispcrsionless [|l5|-j30| and (weakly) dispersive me- 
dia p|,|31 -|36|| . Formulations in which the medium is in- 
troduced macroscopically via a constant (i.e., frequency- 
independent) permittivity may run into problems when 
dispersion is present [ ^5| , ^8p4|j25[ | . The problem of quan- 
tization of radiation in dispersive and absorbing lin- 
ear media has been considered in |37]-|4(|. Recently, the 
method of quantization via Green-function expansion 
ittoj , which is closely related to the (mesoscopic) diagonal- 
ization scheme considered in |39|, has been extended to 
arbitrary linear media, i.e., arbitrarily space-dependent, 

An ex- 
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tension of this approach to nonlinear, absorbing media 
has been given in |l4| ]. An interesting attempt to treat 
the problem of quantization of radiation in dielectric me- 
dia composed of two-level atoms has been considered in 



In this article the propagation of quantum solitons in 
absorbing nonlinear wave guides is studied. It is well 
known that nonlinear waveguides are promising tools 
for "turn-key" generation of nonclassical light Q, be- 
cause of the great nonlinear interaction length and the 
high intensity of the radiation |J. Moreover, the radi- 
ation can be well localized in space and time by form- 
ing solitonlikc pulses. The stability of classical solitons 
in a lossless medium is based on the balance between 
anomalous group velocity dispersion and Kerr-induced 
self-phase modulation |f9|,^0| . The theoretical description 
of quantum solitons usually starts from a Hamiltonian, 
the classical form of which leads to the classical (opti- 
cal) nonlinear Schrodinger equation |p3| , p7| , ^T| , ^2| . The 
quantum nonlinear Schrodinger equation (QNLSE) can 
be solved by means of the Bethe ansatz |3^ , |5l| , |52t . In 
practice, however, the QNLSE must be supplemented 
with further terms in order to describe effects such as 
absorption, third-order dispersion, and Raman scatter- 
ing. Therefore methods based on the Hartree approxi- 
mation []3l|,[53| and linearization around the classical so- 
lution |54|-|58[ have been developed that also apply to the 
solution of modified QNLSE's. 

Alternatively, phase-space methods |39| can be applied 
[Q,[30-62|, starting from the master equation and intro- 
ducing — after appropriate spatial discretization high- 
dimensional phase-space functions that satisfy (pseudo-) 
Fokker-Planck equations (PFPE). The resulting nonlin- 
ear, high-dimensional, partial differential equations are 
hard to solve in general. In the case of proper FPE's the 
solution can be tried to be obtained by stochastic sim- 
ulation. In what follows cumulant expansion (see, e.g., 
[^3|,Q) is used in order to derive an infinite set of evo- 
lution equations for the cumulant hierarchy associated 
with a chosen phase-space function. Introducing appro- 
priate truncation, the resulting finite set of ordinary dif- 
ferential equations can then be solved numerically. The 
method is applied to damped optical quantum solitons in 
absorbing Kerr media, and numerical solutions in Gaus- 
sian approximation are presented. The quantum noise is 
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calculated, with special emphasis on local and spectral 
squeezing fl|,|4|-|lo|| and spectral photon-number correla- 
tion JTJ. 

The article is organized as follows. In Sec. || the under- 
lying master equation is given. The corresponding pseu- 
do-Fokker-Planck equations for s-parametrized phase- 
space functions are derived in Sec. III. Section IV out- 
lines the cumulant expansion and the derivation of the 
cumulant evolution equations. Numerical results of the 
solution in a Gaussian approximation are presented in 
Sec. |\], and a summary and some concluding remarks 
are given in Sec VI. 



II. MASTER EQUATION 

In the simplest model absorption can be described in 
terms of Markovian relaxation theory, introducing in the 
QNLSE a linear damping term and an additive operator 
noise source associated with the losses |^ , |6l| . A more 
rigorous quantization scheme shows that the nonlinear 
interaction can also mix the field and the noise source to 
give a modified QNLSE with multiplicative noise terms 
Jl4| . Only for sufficiently weak absorption can the simple 
model of additive noise be applied. Here we will restrict 
our attention to this case. 

Let us first consider pulse propagation in a lossless 
fiber. When the pulse is long compared to the optical 
period and spectrally well localized in an interval Ak 
around the carrier wave number k c , then the dispersion 
relation uj(k) (of the linear wave equation) can be treated 
in a second-order approximation, 



(1) 



Studying the nonlinear media it may be convenient to de- 
fine the optical (mid)frequency uj c slightly different from 
L0 C given in Eq. (m) 



5lo, Suj <C u) c 



(2) 



Separating from the (positive- and negative-frequency 
parts of the) field the rapidly varying exponentials 
exp[=Fz(d> c i — k c z)], the pulse can be described in terms of 
slowly varying annihilation and creation operators 6(z, t) 
and w(z, t), respectively, that satisfy the equal-time com- 
mutation relations Hi 



b(z,t),&(z',t) = A~ 1 5(z-z'), (3) 



b(z,t),b(z',t) = tf(z,t),fr(z',t) = (4) 



(A is the effective cross section of the fiber) . In what fol- 
lows it will be convenient to use a reference frame that 
moves with the group velocity, i.e., 



(5) 



dk 



(0) 



and to introduce the operators 

a(x,t) — b(v gT t + x,t), a)(x, t) = b)(Vgtt + x, t). (7) 

Obviously, the operators d(x,t) and d^(x,t) also satisfy 
commutation relations of the type given in Eqs. (||) and 
(j|) [with S(x - x') in place of S(z - z')]. 

Under the assumptions made, the undamped spa- 
tiotemporal pulse evolution is governed by the Hamil- 
tonian ip^M^ 



H = hA dx 



5uja)a 



+ \uj^ (g x at) (d x d) + Ix^tfaa 



(8) 



where the constant \ is related to the third-order suscep- 
tibility x (3) as § 



X 



3x (3) ^ g rfcc) 2 

4e?e 



(9) 



(e r is the relative permittivity at frequency ui c ). Note 
that solitons can be formed either in focusing media with 
anomalous dispersion (x< 0, w' 2 ' > 0) or in defocusing 
media with normal dispersion (X>0, <jW<0). 

Including in the theory damping as outlined above, the 
density operator (in the Schrodinger picture) satisfies the 
master equation 



ih d t p = [H, p] + i^Lp, 
where the Lindblad term Ip5l 



(10) 



"fLp = jhA / dx 



JV t h (2d) pa — pad) — da) p) 



+ (^Vth + 1) (2apcr — pa) a — a) dp) 



(11) 



is responsible for damping, with 7 being the damping 
constant, and 



iVth = 



exp 



knT 



- I 



(12) 



(&b is the Boltzmann constant, T the absolute tempera- 
ture). 

To obtain a discrete version of Eq. (|Io|), we divide the 
x interval into m cells of size Ax numbered with index j 
and introduce discrete operators |61[] 



d 3 (t) 




dx a(x, £), 



(13) 



where 



with 
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Xj =Ax(j-[m/2}). (14) 

It follows from Eqs. (||) and (Q) that the operators d.j and 
a\ satisfy the equal-time commutation relations 



aj(t),aUt) 



0. 



(15) 
(16) 



When Ax is sufficiently small, then Eq. ( |l3| ) reduces to 
aj(t) ~ a(xj,t)v'Ax>l, (17) 

and the relations for continuous quantities can be re- 
placed with the relations for discrete quantities according 
to the rules 



dx T\a(x), a/{xf\ 



AxA, a]/vAxA 



d r a 



1 



VAxA Ax 



(18) 



(19) 



where J-[d(x), a*(x)] is an arbitrary function of a(x) and 
d^(x). Finally, it is often convenient to describe pulse 
propagation in terms of appropriately scaled quantities. 
In our numerical calculations we have used the scaled 
quantities given in Appendix [Aj 



III. FOKKER-PLANCK EQUATIONS 

It is well known that the description of the quantum 
state in terms of the density matrix is fully equivalent 



to a description in terms of phase-space functions. In 
what follows we restrict our attention to s-parametrized 
phase-space functions ^,^], which for m operators ctj 
are functions of m complex phase-space variables Oj . For 
notational convenience we will omit the time argument 
and write the phase-space arguments in a compact vector 
form 

P(ai,a 2 , ■ ■ .,a m ,t; s) = P(a; s) = P(u, v;s), (20) 

with the real vectors u and v being defined by 

a = u + iv. (21) 

Note that expectation values of s-ordered operator pro- 
ducts are obtained by averaging the corresponding pro- 
ducts of complex phase-space variables with the s- 
parametrized phase-space function, i.e., 



,trii 



iv m ) 



= J d m ud"W [(ui -ivi) ni ...{v 
x (iti + ivi) kl ...(u m + iv m ) km 1 P(u, v; s), (22) 



where the notation (d\ ni . . . a^") s is used to introduce s 
ordering. 

Applying standard rules (see, e.g., |||), the evolution 
equation for P(u, v; s) can be derived from the master 
equation straightforwardly. After some algebra we obtain 



ftP(u,v; S ) = ^ 

3 

d 



5u> 



dUj 



dVj 



,( 2 ) 



2Ax 2 



dUj 



(2v j -v j+1 -v j - 1 ) 



dv 



■(2uj-u j+ i-Uj-i) 







duj 3 dvj 



d 



d 



duj 3 dvj 



X 2 



dv] 



^3 U 3 



d d 
duj dvj 



? I Nth 



1-s 2 
16 

l-.s 



M + dvi) p(u ' v;s) 



(23) 



[x = xl (AjrA)]. An equation of this type is also called the 
pseudo-Fokker-Planck equation. Apart from the fact that 
in an ordinary FPE derivatives only up to second order 
appear [ |6S| ] , it is impossible, in general, to find a stochas- 
tic process that is equivalent to Eq. ( p3| ) and can be used 
for solving the equation. For s = 1 (P function) and s = 
— 1 (Q function) the two last terms in curly brackets in 
Eq. ( (23[) correspond to noise with a nonpositive diffusion 
matrix which cannot be modeled by a stochastic process 
j|. Nevertheless, such a PFPE can have nonsingular so- 



lutions (see, e.g., J7(|). For s — Q (Wigner function) the 
fifth term in curly brackets in Eq. (^3|) contains so called 
third-order noise which also cannot be treated by stochas- 
tic simulation |7l| . For all other values of s the two terms 
appear together. It should be pointed out that for the 
so-called positive P representation the master equation 
( p0| ) corresponds to a proper FPE with positive definite 
diffusion matrix P,|6l| , [72| , so that stochastic simulation 
applies (see, e.g., ]73[). However, the price to be paid is 
the doubling of variables. 
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IV. CUMULANT EXPANSION 



In order to solve Eq. (p3|), we apply cumulant expan- 
sion techniques. For this purpose we express the charac- 
teristic function of P(u, v; s) in terms of cumulants and 
derive the cumulant evolution equations. After appropri- 
ate truncation the set of ordinary differential equations 
can be solved numerically. 



A. Basic definitions 

Let us consider a (quasi) probability distribution func- 
tion P(q) = P(qi, Q2, ■ ■ ■ ,qf) of / real variables qu and 
define its characteristic function x(q) =x(q,, q 2 , • ■ ■ , q ) 
by the Fourier transform of -P(q), 



P(q) = J^jf J d/ Sx(q)^- q , 



(24) 



x exp 



W ■■■q} 1 )) = 



x In 



{M 







■■If 


»£>: 


..k f l 


-) 


Til 












■ 


■iA 


hi. 


•A/1 



q=0 



q=0 



(30) 



(31) 



{{{q® ■ ■ ■ <$)) = 0). Note that the symmetry relations 



rii «j 



rij rii 



(32) 



x(q) = / df q p( q y^. 



(25) 



From Eq. ( |25| ) it follows that the moments (q™ 1 . . 

of the distribution function P(q) can be obtained from 

x(q) by differentiation, 



(q^...q n /)= J dfqq?...q n /P(q) 



idq l 



x(q) 



(26) 



q=0 



(rifc > 0). Hence, the Taylor-series expansion of the func- 
tion x(q)i which can be regarded as moment generating 
function, reads as 

(iq ) ni ...(iqTf 

x(q)=E<tf i -g/ / > 1, I , • ( 2? ) 

— z — ' 1 Til . . . Ufl 

{n k } 1 

Introducing the cumulant generating function $(q), 

x(q) = e*<3>, (28) 

the cumulants ((g™ 1 . . . q^)) are defined by the Taylor- 
series expansion of $(q) as follows: 

(io ) ni (iq ) n f 

$(q) = 2 iii 1 ■ ■ ■ <z?» _ ,"_T ■ ( 29 ) 



m! . . . n./! 



Combining Eqs. (27) - (|29|), the cumulants and moments 
can be related to each other as 



(33) 



are valid. For a Gaussian distribution higher than second- 
order cumulants are equal to zero, 



(Qi 1 



If 



if 



n f >2, 



(34) 



so that higher than second-order moments can be ex- 
pressed in terms of first- and second-order cumulants 
((q,)) and ((ft£fr», respectively (see, e.g., [nj]). 

Obviously, Eqs. (24) - ( |33| ) also apply to the s-para- 
metrized phase-space functions P(u, v; s) introduced in 
Sec. @ [P(q) P(u,v;s), x(q) -> x(u,v;s), $(q) -> 
<5(u, v;s)]. Since the characteristic functions for differ- 
ent values of s are related to each other as |67) 

X(u,X;s) = x(u,v; s') 

x exp[i(s - s ')(ul + vl + ... + ul + , (35) 



from Eqs. ( J28| ) and (|29|) it is seen that only the second- 
order cumulants of equal variables depend on s, 



((UjU r )) s = ({UjUj>)) a > - Us - s')5j 



(36) 



and 



accordingly. All other cumulants are inde- 



pendent of s [|63| . Note that the mean value of the com 
pie 
by 



plex amplitude (a,j) and the intensity Ij = (ctja,-) are given 



(%) = ((Uj))s+i((Vj))s, 



u J//s 



(37) 

+ Us-1). (38) 
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B. Evolution equations 

Let us assume that a (quasi)probability distribution 
function P(q) satisfies an evolution equation of the fol- 
lowing type: 



5 t P(q) 



= E 

{n,},{h} 

d xni 



fei 



Ak x ...k f 

JlJll. . .71 f 



d 

dq f 



k f 
If 



P(<l), 



(39) 



where An^'.'.Jj are constants. It follows from Eq. j39| ) 
that the characteristic function x(q) ; Eq. (|25|), satisfies 
the evolution equation 



dtX($ = d t J d / qe^ q P(q) 

.ki.-.kf (_■ \m 
-fi-n 1 ...n f [ "f^J 

\ k i / \ fc f 

\ I d \ 



E 

{ni},{ki} 



x V d lJ 

and hence [see Eq. (p8|)] 



x(q), 



(40) 



5t*(q) 



-*(q) 

{nj.ffcj 

/ 9 \ 



4 1 1 ;:.n / / (-%) ni ---(-i£ / )" / 

1 ( d \ k , 



idq 1 



idq 



e *(q) 



(41) 



Substituting for $(q) in Eq. ( |4l| ) the Taylor-series expan- 
sion (29) and comparing the coefficients of equal powers 
we arrive at 



A k i--- k f 



q r /))= 

{ni},{h} 

(_j)™i+fci+»'i ri i (—i) n f +k f +r frf\ 



*(q) 




(42) 



q=0 



(l/rt! = if n<0). Performing in Eq. (|42j) the differentia- 
tions yields the cumulant evolution equations explicitly. 

We now apply the scheme to the PFPE ( p3| ) in order 
to obtain the evolution equations of the cumulants asso- 
ciated with P(u, v; s). For the sake of clearness we some- 
what change the notation, using C„„... „t.. v (i, x' , . . .) to 
denote the cumulants (for chosen s), e.g., 



i (x, x') = ({u j v j u],)) a (j ^ /), 



(43) 



with x = Xj [cf. Eq. (|14|); note that Xj±% —x±Ax]. After 
some algebra the evolution equations for the first-order 
cumulants C u (x) and C v (x) are derived to be 

d t C u (x) = —jC u (x) +SujC v (x) 

+ 7T7-^ [2C„(x) - C^x-Ax) - C„(:r+Ax) 
2Air z 

+ xC v (x) [Cl{x) +C 2 v (x)] 

+ xC v (x) [{a - 1) + C„ 2 (x) + 3C v 2(x)} 

+ 2 x C u {x) C uv (x) + x [C u 2 v (x) + C v3 (x)] , (44) 



dtC v (x) = — 7 C„ (x) - <Sw C„ (x) 

f9l 



- — -=■ [-2 C„(x) + C„(;r-Ax + C u (x+Ax)] 

-xCu(x) [Cl{x) +C*(x)] 
- X C u (x) [(s - 1) + 3C U 2 (x) + C v2 (x)] 
2 X {x) C uv (x) - X [C U 3 fx) + C UV 2 {x)\ . (45) 

Owing to the nonlinear interaction the first-order cu- 
mulants are coupled to second- and third-order cumu- 
lants, whose evolution equations are rather lengthy [for 
the evolution equations for the second-order cumulants, 



see Eqs. ([bT]) - (Q in Appendix B . N ote that the right- 
hand sides of Eqs. (0), (|45[), and (p32|) are s independent 
[in accordance with Eq. (|36|) and the statement given 
there] . Note that in the limiting case Ax — ► the cumu- 
lant evolution equations [cf. Eqs. (^i|), (pij), and (Bl)- 
( |B3|) 1 are transformed into partial differential equations 
for corresponding cumulant correlation functions. 



C. Gaussian approximation 

With regard to a numerical treatment, the hierarchy of 
cumulant evolution equations, which represents an infi- 
nite set of coupled, nonlinear, ordinary, first-order differ- 
ential equations, must be decomposed into finite, closed 
systems of equations. For this purpose approximation 
schemes are required to be applied. Note that only for 
the linear system (x = 0) does the infinite set of equa- 
tions decompose into finite, closed subsets, each of which 
contains cumulants of the same order. 

For a sufficiently high peak p hot on number of the order 
of magnitude of 10 9 [Eq. (A9)] an expansion in 1/n 
may be performed, in which the second-order cumulants 
are considered to be of magnitude 1/n in comparison with 
the first-order cumulants and cumulants of higher than 
second order are neglected. Comparing terms of equal or- 
der we obtain a linear system of differential equations for 
the second-order cumulants, in which the first-order cu- 
mulants governed by truncated Eqs. (E5|) and ( f45| ) play 
the role of time-dependent coefficients (for such lineariza- 
tion schemes for lossless systems, see ||,[34-57 , and for 
absorbing case, e.g., fl58|). The approximation scheme 
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may be extended straightforwardly in order to also in- 
clude higher-order cumulants in the calculation. 

Another approach to the problem is to retain all cu- 
mulants up to some appropriately chosen order Af and 
to set the remaining higher-order cumulants equal to 
zero. In what follows we restrict our attention to the 
second-order (Gaussia n) a ppr oxim ation (TV = 2). Thus 
in Eqs. @, (§f), and (jB|) - (Hf) all terms that contain 
cumulants of third and fourth order are omitted, and the 
now closed system of nonlinear differential equations is 
solved numerically. This approximation obviously corre- 
sponds to the assumption that the quantum state of the 
pulse can be approximated by a Gaussian phase-space 
function P(u, v; s). Note that for Gaussian distributions 
higher than second order cumulants exactly vanish. 

The Gaussian approximation may be expected to ap- 
ply if the damping is not too small. To roughly estimate 
the range of validity of the Gaussian approximation, we 
note that the linearization approach used in J54|] for an 
undamped fundamental quantum soliton fails at (scaled) 
times t^n 1 / 4 |52| (for scaled quantities, see Appendix [A]) . 
On the other hand, in absorbing media the soliton noise 
relaxes to thermal-equilibrium noise on a time scale ~ 
7 _1 . As the non-Gaussian correlations are suppressed 
more rapidly than the Gaussian ones it is expected that 
for 7" 1 < n 1 / 4 the Gaussian approximation is valid for 
all times of pulse propagation. Hence, the lower limit for 
the scaled absorption coefficient is given by 



7> 



-1/4 



(46) 



which for n~ 10 9 yields 7>5.6 x 10 -3 . For smaller values 
of 7 the Gaussian approximation is expected to be valid 
only in a limited time interval t^n 1 / 4 . 



V. RESULTS 
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In the Gaussian approxi matio n t he s teady-state solu- 
tion of Eqs. (|J), (11), and @) - @ can be found to 
be 

C u (x) = C v {x) =0, (47) 
C UtU {x, x') =C ViV (x,x') = ±[N th + ±(l- a)] S(x, x'), (48) 
C UtV {x,x') = C VtU (x,x') = (49) 

[S(x, x') = l, i{\x — x'\< Ax/2, and S(x, x') = otherwise]. 
It corresponds to thermal fluctuations, without any cor- 
relation between the fields at different space points. In 
the numerical solution of the time-dependent problem we 
have used the steady-state values ( |48|) and (fl9|) as initial 
values of the second-order cumulants. With regard to the 
initial values of the first-order cumulants, we have used 
the fundamental soliton solution of the classical nonlinear 
Schrodinger equation p9|,p0|, 



C u (x) = y/n^ sech(a;/xo), 
C v (x)=0 



(50) 
(51) 



FIG. 1. The spatiotemporal evolution of the intensity 
/ of a soliton is shown for various values of the damping 
constant 7. 



[no = uAx/xq, Eq. ( |A8| )]. The initial conditions are re- 
alized by a multimode displaced thermal state, without 
entanglement between the modes corresponding to differ- 
ent space points (for displaced thermal states, see, e.g., 

To be more specific, we have performed the calcula- 
tions assuming a peak photon number of the initial pulse 
of n = 10 9 and a reservoir photon number of N t h = 10~ 16 
and setting the parameter 5u equal to zero. The value 
of Nth corresponds to a (vacuum) carrier wavelength of 
A c = 2ttc/uj c — 1.5 /im of the pulse in vacuum and a 
temperature of T = 300K [see Eq. @]. Thus the pulse 
is initially prepared in a displaced thermal state that is 
almost a coherent state. Assuming, e.g., losses of T = 



G 



0.3dBkm _1 [Eq. (A6)] and a fib er dispersive parameter 
of D= 20psnm~ 1 km -1 [Eq. JA5| )1, the (scaled) damping 
constant is 7 = 5.8 x 10 -3 for a pulse duration of to = 
2 ps, and 7 = 1.4 X 10 _1 for tg = 10 ps, and the dispersion 
lengths are Xd = 170 m and Xd = 4.2 km, respectively [see 
Eqs. @ - ©]. 




FIG. 2. The local field fluctuations of an undamped 
soliton (7 = 0) are shown for fixed propagation time t 
and various positions x in the pulse. The uncertainty 
ellipses in the phase space (s = 0.85) are plotted in ar- 
bitrary units around the (solid line) complex amplitude 
{a{x)) = {u{x))+i{v{x)). 



We have solved the equations numerically on a grid of 
to = 200 points with Arr/rro = 0.1 and absorbing bound- 
ary conditions, employing the Runge-Kutta method of 
order 8(5,3) with step-size control []75[. In order to test 
the program, we have performed the calculations for dif- 
ferent values of s and verified Eq. ( p6| ) to be within the 
numerical round-off error. Examples of the temporal evo- 
lution of the pulse intensity / = I(x) = Ij, Eq. (|3^), are 
shown in Fig. ^ for various values of 7. 



with Eq. ( |36[ ) it can be seen that b and B depend on s, 

b^b s = b Q -\s (54) 

(and B accordingly), whereas ip is s independent. The 
positions in phase space of the ellipses are determined 
by the s-independent first-order cumulants C u {x) and 
C v (x). Squeezed local fluctuations are observed if b is 
reduced below the vacuum level & vac , 



b < b v 



(55) 



(Appendix |^). Obviously, the property of a state to be 
squeezed does not depend on s. For illustrational reasons 
it is convenient to choose s as large as possible [otherwise 
the difference b— (1 — s)/4 may not become visible]. 



t = 8i 



t = n 



t - 6t 
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FIG. 3. The temporal evolution of the local field 
fluctuations of an undamped soliton (7 = 0) is shown 
for three typical positions in the pulse. The uncertainty 
ellipses in the phase space (s = 0.85) are plotted in ar- 
bitrary units around the complex amplitude (&(x)) = 
(u(x)} + i{v(x)). For comparison, the vacuum noise is 
indicated by circles. 



A. Local fluctuations 

In the Gaussian approximation, the local field noise 
(expressed in terms of s-ordered variances) can be visu- 
alized in phase space by means of uncertainty ellipses 
with the large B 1 / 2 and small 6 1 / 2 axes and orientation 
angle ip, 



B.b 



C U 2+C V 2± J(C u 2-C v i) Z +4C^ 



2cp = arg (C u 2 - C v 2 + 2iC uv ) 



(52) 



(53) 



(Appendix |C|) , where the argument x is skipped for sim- 
plicity. From inspection of Eqs. (151) and (m) together 



The spatiotemporal evolution of the local field fluctu- 
ations (together with the complex amplitude) of an un- 
damped soliton (7 = 0, x > 0' w ' 2 ' < 0) i s illustrated 
in Figs. ^ and |[ Since the uncertainty ellipses are very 
small (^^o 1 ^ 2 ) compared with the (absolute) values of 
the complex amplitude (~1), the ellipses are plotted on 
an enlarged scale (of arbitrary units). Figure || represents 
the dependence on space of the local fluctuations for fixed 
time. The vacuum-noise level can be estimated from the 
uncertainties observed sufficiently far from the pulse cen- 
ter (| a; I ^ 4xo in the figure), where the field is prepared 
in a state close to vacuum. The Kerr-induced self-phase 
modulation may be regarded as being the main reason 
for the enhanced phase noise of the field in the vicinity 



7 



of the pulse center (for continuous wave radiation, see, 
e.g., |7q| ). The reason for the spatial redistribution of 
the local fluctuations may be seen in the dispersion that 
also leads to the observed enhanced amplitude noise at 
the wings of the pulse (for the influence of the dispersion 
on the squeezed pulses in linear media see, e.g., |77]]). In 
Fig. H the temporal evolution of the local fluctuations is 
shown for various space points. It is seen that during the 
pulse propagation the complex amplitude rotates clock- 
wise (counter-clockwise in the case \ < 0; > 0) m the 
phase space due to the nonlinear phase shift. The period 
of the rotation is close to the classical period 47rid , 
because of the large value of n. For small propagation 
times the field in the pulse center is squeezed, whereas at 
the wings the noise exceeds the vacuum level (indicated 
by the circles in the figure). In the further course of time 
an increase of the noise is observed. 
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FIG. 4. The temporal evolution of the minimum-qua- 
drature noise b in the center of a soliton is shown for 
various values of the damping constant 7 (s = 0). 



The influence of absorption on the spatiotemporal evo- 
lution of the minimum quadrature noise ~ 6 is shown in 
Figs. |J and ||. From Fig. Q that presents the temporal evo- 
lution of 6 in the pulse center it is seen that absorption 
reduces the maximally achievable amount of squeezing. 
Figure || reveals that for small propagation times (i < id) 
squeezing is observed over the whole pulse, whereas in 
the further course of propagation (i >id) "noisy" wings 
appear owing to dispersion-assisted noise redistribution 
in the pulse. It is also seen that absorption reduces the 
effect of noise enhancement at the wings. 

The local field noise can be described in terms of 
squeezed thermal states, and the squeezing condition (pi 
can be rewritten as 



[Eq. (|C9|)] is the associated squeeze parameter. The pa- 
rameter r is a measure of the "deformation" of the un- 
certainty ellipses, whereas n determines their "size" . Ex- 
amples of the spatiotemporal evolution of n and r are 
shown in Figs. |^ and |7[ respectively. It is seen that for 
small propagation times r increases more rapidly than n 
in the pulse center, and hence the field in the pulse cen- 
ter can be squeezed after the pulse has entered the fiber. 
The destructive influence of absorption is demonstrated 
in Figs. |(6) and |(6). 



7 = 



lo gio[ fo / & v 




log 10 [6/6 v 
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FIG. 5. The influence of absorption on the spatio- 
temporal evolution of the minimum-quadrature noise b 
of a soliton is illustrated for three values of the damping 
constant 7 (s = 0) . 



(2n + l)e 



-2r 



< 1 



(56) 



(Appendix |(|). Here, n [Eq. (C8)] is the "thermal" pho- 
ton number associated with the local field noise, and r 
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°-4 



0-4 



7* d = 0.05 




FIG. 6. The spatiotemporal evolution of the "ther- 
mal" photon number n of the squeezed thermal state 
associated with the local field noise of a soliton is shown 
for two values of the damping constant 7. 



jt d = 0.05 




FIG. 7. The spatiotemporal evolution of the squeeze 
parameter r of the squeezed thermal state associated 
with the local field noise of a soliton is shown for two 
values of the damping constant 7. 



B. Squeezing spectrum 

The squeezing spectrum of the (output) pulse can be 
measured using balanced homodyne detection J62]. Ne- 
glecting outcoupling effects, the temporal profile of the 
pulse after propagation through a fiber of length L is 
described by the operator 



a(L — v gr t, t) — a{— rugr, t + L/vgf), 



(57) 



where r = t — L/v gr is the difference of the absolute time 
t and the pulse arrival time L/v gT . The approximation in 
Eq. (Q) requires the temporal pulse width t n to be small 
compared with the dispersion time id, Eq. (A3), 



to/t d ~ Afc|w< 2 >|/|u; (1) l «1. 



(58) 



Hence, the dependence on r of the second argument in 
Eq. ( p7j ) can be ignored, and the time profile of the pulse 
can be described by 



o(— TV gI , L/Vgr) — a(x, L/v gI ) = d(x), 



(59) 



where for notational reasons the time coordinate r is for- 
mally replaced by the spatial coordinate x 

From Appendix [j5], the squeezing spectrum S(oj) can 
be expressed in terms of (discrete) s-ordered cumulants 
as 

S(u) = 2Re[f L (- W , w ) + e 2 ^G L (- W , w )] Iq\ (60) 
where 

t-L\u, u ) - — \a Lj a Lj >e 1 1 
33' 

X {((u jUj >)) s + ((VjV f )) s + |(s - 1)%- 

-i [{{ujVj>)) a - ({vjUj>)) a ] } , (61) 



ii' 



X {{(ujUf)) s - {{VjVy)), 



} 



(62) 
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Here, Jo is a normalization constant and o,l j is the com- 
plex local-oscillator (discrete) amplitude, with (p being 
the relative phase between the local oscillator and the 
signal. Negative values of S(cu) signify squeezing at fre- 
quency u>. Note that the squeezing spectrum involves in- 
ternal pulse correlations expressed in terms of the two- 
argument cumulants in Eqs. ( |6l| ) and ((6^). For the chosen 
frequency lj the minimum of S(lo), 

S min (uj) = 2[ReF L (-uj,Lo) - |G L (-w, u)\] If 1 , (63) 

is observed if the phase tp is chosen such that 

2^ = 7r-arg[G L (-w,w)]. (64) 



t/t d 




FIG. 8. The influence of absorption on the temporal 
evolution of the midcomponent of the squeezing spec- 
trum S m in(0), Eq. (ppj), is shown. 



The temporal evolution of the minimum squeezing 
spectrum S m i n (uj) is shown in Figs. |§| and^ (ljq = 2/xo). 
It is assumed that the local oscillator (LO) pulse has the 
shape of the initial fundamental-soliton pulse. Figure |§| 
demonstrates the influence of absorption on the temporal 
evolution of the squeezing spectrum for lo — 0, and exam- 
ples of the full squeezing spectrum are shown in Fig. |^. 
At the initial stage of propagation wideband squeezing 
is observed which evolves to enhanced sideband noise in 
the further course of time [Figs. ||(a) and 0(6)]. From 
Fig. ||(a) it is seen that in the lossless case the sidebands 
increase with time, whereas the midcomponent decreases 
(see also Fig. . Absorption introduces bounds, as can be 
seen from Fig. [)](&). In particular, for sufficiently strong 
absorption the sidebands do not appear at all [Fig. |9](c)]. 

It should be noted that in the numerical implementa- 
tion of the x integrals the value of Ax imposes according 
to the sampling theorem [7c] ] a condition on the upper 
bound of frequency: \u)\ <cJ max =7r(Aa:)~ . Further, there 
is a minimally resolvable frequency difference Aw m i„ = 
2ir(mAxy 1 . For Ax/x Q = 0.l and m = 200 wc find that 
Wmax = 57T oj and A = 0.057ro;o. 



C. Spectral photon- number correlations 

Recently experiments have been performed in order to 
measure the correlation of the number of photons N(lo) 
at different frequencies [[nj. A measure of the degree of 
correlation is the correlation coefficient 




jt d = 0.05 




FIG. 9. The temporal evolution of the squeezing 
spectrum S m i n (us), Eq. (^), is shown for various values 
of the damping parameter 7. 
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7 = >yt A = °- 05 





FIG. 10. The photon-number correlation coefficient 
T](u>,u}'), Eq. d&5|) , of an undamped soliton (7 = 0) is 
shown for two values of the propagation time t. Note 
that the symmetry relation r)(u),u') = 77(0/, w) is valid. 



FIG. 11. The photon-number correlation coefficient 
r](u>,u'), Eq. (|65|), of a damped soliton (7^ = 0.05) is 
shown for two values of the propagation time t. Note 
that the symmetry relation r](u),ui') = rj(u)',ui) is valid. 



where (: AN(ui)AN(ui') :) is the normally ordered pho- 
ton-number covariance, which is related to the ordinary 
covariance as 

(AN(uj)AN(u;')) 
= (: AN(u))AN(u)') :) + (N(u)))5(w, uj') (66) 

(Appendix |e|). The second term in Eq. (|66|), which rep- 
resents a (quantum) shot-noise contribution, vanishes for 
lu =/= lu' , and hence t](lu,lu') is the ordinary correlation 



coefficient. For lu = 



the autocorrelation coefficient 



t)(uj,lv) is nothing but the (normalized) spectral photon- 
number variance in normal order. Negative (positive) val- 



ues of f](u>, to ) correspond to sub-Poisson (super-Poisson) 
photon- number statistics at chosen frequency lu. 

Expressing (N(us)) and (: AN(ui) AJV(w') :) in terms of 
(discrete) s-ordered cumulants yields 

(N(oj)) = \f(-w,u) + \E(uj)\ 2 ] , (67) 

(:AN(lu)AN(lu'):) = \F(-u, oj')f + \G(lu,lu')\ 2 
+ 2Rc [F(-w, lu')E(lu)E*{lu') + G(w, w')E*{u>)E*{u)')] , 

(68) 

where 
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FIG. 12. Contour plot of the photon-number cor- 
relation coefficient tj(lj,lj'), Eq. (pa), of an undamped 
soliton for the same parameters as in Fig. [l^. 



FIG. 13. Contour plot of t he p hoton-number corre- 
lation coefficient rj(u,LJ r ), Eq. (p5|), of a damped soliton 
for the same parameters as in Fig. [llj 



'AxAo 
2tt 



E 



Ul )) s +i v 



3 lis i *\\"3//sl ) 



jj' 



(69) 



+ ({ujUj')) s + {{vjV f )) s +i [((%^-')} s -((w)) s ] }, (70) 



- + i [({ujVj>))s + {{vjUj>)) s ] }, (71) 

with Aw being the resolving frequency (see Appendix ||]) . 
We have performed the calculations choosing the mini- 



mally possible value of Au in the d iscre tization scheme 
used (see the last paragraph in Sec. VB ). 

In Figs. [l(] and [jj] the correlation coefficient rj(uj, a/) of 
an undamped soliton is plotted for two values of the prop- 
agation time t. The figures show typical features mea- 
sured recently [fll|| . In particular, in the central part of 
the spectrum negative correlations are observed, whereas 
outside the center the frequency components can be 
strongly positively correlated. The negative values of the 
autocorrelation coefficient r)((£, oj) (solid lines in Fig. |l(i| ) 
reveal, in agreement with Q, spectral sub-Poissonian 
statistics in the central part of the spectrum followed on 
both sides by regions of super-Poissonian statistics (dot- 
ted lines in Fig. [10| indicate the shot-noise level). With in- 
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creasing propagation time the cross-correlation strength 
and the super-Poissonian autocorrelation strength are in- 
creased, whereas the sub-Poissonian autocorrelation ef- 
fect becomes weaker [compare Figs. |l(](a) and |l2](a) for 
t — 2.5 id with Figs. 10(6) and 12(6) for t = 5 id respec- 
tively] . 

In Figs, [ll] and [l^ the correlation coefficient rj(co,u)') 
of a damped soliton is shown for the same values of the 
propagation time as in Figs. [l(] and [l2| Whereas in the 
early stage of propagation similar results are observed 
[compare Figs. |l(](a) and |l2|(a) with Figs. [n](a) and 
|l3|(a), respectively], in the further course of time both 
the negative and positive cross correlations are smeared 
owing to damping [compare Figs. |h](6) and [l2](6) with 
Figs. |ll](6) and |l^(6), respectively]. In contrast to the 
cross correlations, the effect of sub-Poissonian autocorre- 
lations in the center of the spectrum can be stronger for 
the damped soliton than for the undamped one. This sur- 
prising effect may be a result of damping out destructive 
fluctuation interferences that occur for vanishing damp- 
ing. 

As mentioned above, from the theory sub-Poissonian 
statistics is expected to occur in the central part of the 
spectrum. Further the calculation yields a correlation co- 
efficient that fulfills the symmetry relation r)(u>,u>') = 
r](—uj,—u)'). The two effects have not been observed in 
the experiment in p| . A reason for the discrepancies 
may be seen in phonon-scattering-induced noise which is 
not considered here. It should be mentioned that for small 
damping an observable sub-Poissonian effect is expected 
to occur only at a rather early stage of pulse propagation 
[Figs. M{a) and 0(a)]. 



VI. SUMMARY AND CONCLUDING REMARKS 

We have studied the influence of absorption on the 
quantum statistics of optical solitons in Kerr media, in- 
troducing multivariable s-parametrized phase-space dis- 
tributions for describing the multimode quantum state 
of a solitonlike pulse. The evolution equations for the 
phase-space distributions can be tried to be solved us- 
ing cumulant-expansion techniques. When the damping 
is not too small, then the quantum state can be treated 
in Gaussian approximation; i.e., only cumulants up to 
second order are to be taken into account. Otherwise, 
the evolution time must not be too large. Starting with 
an undamped fundamental soliton prepared in a nearly 



coherent state, we have solved the set of differential equa- 
tions in Gaussian approximation numerically. From the 
solution, we have analyzed the local field fluctuations in 
terms of minimum-quadrature noise and squeezed ther- 
mal states, and we have calculated the squeezing spec- 
trum and the spectral photon-number correlation coeffi- 
cient. 

The cumulant evolution equations in the rigorous 
Gaussian approximation used in this article contain non- 
linear terms that are disregarded in the linearization ap- 
proximation. The obtained results show that for suffi- 
ciently big coherent amplitude the Gaussian approxima- 
tion closely corresponds to the linearization approxima- 
tion. Deviations of the exact behavior of the quantum 
noise from the behavior obtained in the linearization ap- 
proximation may therefore be regarded as being of non- 
Gaussian type. In this way, the cumulant method offers a 
possibility of improving the theory, including in it higher- 
order cumulants — at least third-order cumulants — for 
describing the non-Gaussian effects in a consistent way. 

The method can of course be used in order to study 
further aspects of damped quantum soliton motion. For 
example, the class of initial conditions can be extended. 
In practice, the pulses which are fed into a fiber are not 
of the fundamental-soliton type in general. Further, se- 
quences of pulses can give rise to interactions between 
them. When the losses are very small, then non-Gaussian 
effects may be observable, which requires inclusion in the 
theory of higher than second-order cumulants. The model 
of the Kerr medium can also be improved by including ef- 
fects like Raman and Brillouin scattering |3lJ , third-order 
dispersion |53|j , and frequency-dependent absorption. In 
order to further specify the nonclassical features of soli- 
tonlike pulses under realistic conditions and compare the 
quantum motion with the classical one, additional prop- 
erties must be studied. With regard to the multimode 
structure of solitons, the study of nonclassical internal 
pulse correlations plays an important role. Finally, appli- 
cation to other systems is also possible, e.g., solitons in 
media |6^], self-induced transparency solitons Jr9| , 
or video pulses M . 
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APPENDIX A: SCALED VARIABLES 



The scaled coordinates x and t are defined by 



x 

x = — 

x 

t 

t = T d - 



(Al) 
(A2) 
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where xq is the spatial pulse width, and t<j is the disper- 
sion time, which is related to the dispersion length Xd 
as 



Xd 

£d = — 



and Xd can be given by [bOl 



Xd 



2 

x v gr 



|/fc( 2 )| \lu^\ 



where to = xq /u gr is the temporal pulse width and 



*«0 = -^ = -2.fS 



(A3) 



(A4) 



(A5) 



[for w( 2 ) and u gr , see Eqs. (|T|) and (§)]. In Eq. ©, the 
(fiber) dispersive parameter Z? is introduced |50(| . Further 
the scaled damping parameter 7 is defined by 

7 = 7 t d = 0.05 • lnlO ■ T [dBkm -1 ] ■ x d [km]. (A6) 

Finally the scaled amplitude operators 



are introduced, where 



no = n ■ 



Ax 

Xq 



Here 



/ 2 U 



(A7) 



(A8) 



(A9) 



is a measure of the initial photon number of the pulse 
P,|6l| and is typically of the order of magnitude of 10 8 - 
lO^or solitons in optical fibers M. 



APPENDIX B: EVOLUTION EQUATIONS FOR 
THE SECOND-ORDER CUMULANTS 



Using the notation introduced in Eq. (j43|), the evolu- 
tion equations for the second-order cumulants read as 



d t C Uia (x,x') = [N th + 



1- a 



+ s\ [C u (x) C v (x) + C uv (x)] > S(x, x' 



-2~fC u ^ u (x,x') + 5wC u>v (x,x') + 5ljC ViU (x,x') 



CO 



(2) 



2Ax 2 
2 Ax 2 



[C UiV (x,x' - Ax) - 2C u , v (x,x') + C u>v (x,x' + Ax)] 
[C v>u (x - Ax,x') - 2C v>u (x,x') + C ViU (x + Ax,x')] 



+x{c u , v (x,x') [Cl(x')+3CZ(x')]+C v , u (x,x') [C 2 u {x)+iCl{x)\ 

+C u<v {x,x') [(a - 1) + C u 2(x') +3C v 2(x')} + C ViU (x,x') [(a - 1) + C u 2(x) +3C„ 2 (i)] 
+2 G u<u (x, x ) [C u (x) C v (x) + C u (x ) C v (x') + C uv (x) 4- C uv (x )] 
-\-C v (x) C u 2 u (x,x ) + C v (x ) C u<u 2 (a;, x ) + 2 C u (x ) C u ^ uv (x, x ) 
+3 C v (x ) C u v 2 (a;, x ) + 2 C u (a;) C U v,u(x, x ) + 3 C^a;) C„2 )U (x, a; ) 

+Cu,M 2 u( a; i 2; ) + C u v 3 [x, x ) + C u 2 v u (x, x ) + C v 3 >u (x, x ) j> , (Bl) 

dtC UjV (x,x') = x [Cu 2 ( x ) ~ Cv 2 ( x )] S(x,x') - 2jC u , v (x,x') + 6uC v , v (x,x') - Suj C u>u (x, x') 
+ \x [C 2 v (x)+C v 2(x)-Cl(x) - C u 2(x)] 5(x,x') 



;(2) 



2Aai 2 



,(2) 



2Aa; 2 



) +C„ )U (a:,a; / + Aa:)] 
') +C v>v (x + Ax,x')] 



+x{-CuAx,x') [3Ct(x')+Ct(x')] +C v , v (x,x') [C 2 u (x) + 3C 2 v (x)] 

-C UtU (x,x') [(a- 1) + 3C* M 2(a;') + C^(x')] + C VtV (x,a/) [(a - 1) + C u 2 (x) +3C v2 (x)} 
+2 (7^(3;, a;') [C u (a:) C v (x) + C ul , (a;) — C u (a/) C v {x') — C uv (x')\ 
—C u (x) C u v 2 (x, x ) + 2 C u (x) C UVtV (x, x') +3 Ct,(a:) C v 2 v {x, x') 

u 3 (a;, a;') — C uuv 2 (x, x ) + C u 2 vv {x, x') + (7^3 „ (a;, a: ) } , (B2) 
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d t C v<v {x,x')= 1 7 ^JV th +_ — £j - flx[C tt (a?)C fl (a?) + C uv (x)}\ 5(x,x') 
—2^C ViV (x,x') - 5ujC u>v (x,x') - 8ujC v>u (x,x') 
+ 2A^ l C "A x ~ Ax > x ) ~ 2 C uA x , x ) + C u,v( x + Ax, a:')] 

W (2) 

+ r-r 9 [C v , u (a;,a; - Ax) - 2C VtU (x,x) + C VtU (x,x + Ax)} 

- X {C U , V { X , x ') [3 C 2 u {x) + C 2 v (x)] + C v , u (x, x 1 ) [3 C 2 u (x') + C 2 v {x')] 

+C a , v (x, x') {(s - 1) + 3 C u i (x) + C V 2 (x) } + C VtU {x, x 1 ) [{s - 1) + 3 C u 2 (x') + C v i (x')\ 
+2 C v ^ v {x, x ) [C u (a?) Ct, (a;) + C u {x ) C v {x ) -+- C uv (x) + C uv {x )] 
+2 C v (x) C uv ^ v {x, x ) + 2 CV, (a; ) C VtUV (x, x ) + C u (x) (7^2 jt , (x, a; ) 
+3 C u (x) C u 2 „ (a;, a;') + 3C„(a;') C V)U a (a;, x') + C u (x') C VtV 2(x,x') 

+C« 3 ,u( a; 7 x ') + Cu,u 3 (^i 2;') + C v 2 u v {x : x') + C v v 2 u (x, a;')} . (B3) 



APPENDIX C: SQUEEZED THERMAL NOISE 



As n must be non-negative, from Eq. flC8|) it follows that 



For chosen x the (single-mode) s-parametrized charac- 
teristic function x(u, u; s) of the local noise in Gaussian 
approximation, 

x(m,^; s) = ex p[~5 (C U 2U 2 + 2C uv uv + C V 2V 2 )] , (CI) 
can be rewritten as 

X(u,v;s) = exp[-i (By! 2 + by! 2 )] , (C2) 

u' + «;' = (u + i«)e < *\ (C3) 

where B, b, and <p, respectively, are given in Eqs. ( p2| ) and 
(|53|). The characteristic function x(w, s) corresponds to 
a squeezed thermal state 

P = S(Of>t h (n)SHO, (C4) 
which is obtained by applying the squeeze operator 

£(0=exp[i(ra 2 -^ t2 )] (C5) 

on a thermal state pth(n) of mean photon number n (for 
squeezed thermal states, see, e.g., |74|). The second-order 
cumulant parameters in Eq. (CI) can be show n to be re- 
lated to the parameters £ = 7 - e 75 and n in Eq. (C4) as 



C U 2 + C V 2 = (n + |) cosh(2r) - \; 



(C6) 



2iC u 



e i9 (n+i)sinh(2r). (C7) 



Inverting Eqs. (C6) and (|CTj) , the parameters r, 9, and n 
can be expressed in terms of the second-order cumulants 
C u 2, C uv , and C v 2 as, on recalling Eq. (p3), 



B + s/4 



= arg (C„2 - C u 2 - 2iC uv ) . 



(C8) 

(C9) 
(CIO) 



y/(B+\s) (b+is)>{ 



(Cll) 



which is nothing but the Heisenberg uncertainty principle 
[note that according to Eqs. (36) and ( (52| ) the quantities 
B + s/A and b + s/A do not depend on s]. In the squeezing 
condition ( |55| ) for the local field noise the vacuum level 
6 vac can be obtained from Eqs. (C8), with n = and B = 
b = b vac . From Eqs. ([!§) and flC7|) together with Eo, (||) 
it is easy to show that the squeezing condition ( p5[) can 
be rewritten to obtain the condition (^6|) in terms of n 
and r. 



APPENDIX D: DERIVATION OF EQ. (60) 



Following []62j, we introduce the generalized detection 
operator 



X(x) = e iv c(x) + e- iv cJ(x) 
and its Fourier transform 



(Dl) 



where 



X(u) = {2n)- 1/2 J dxX(x)e lwx 
= e^c(w) + e^c t H, 

c(x) = a t (x)a L (x) - (a 1 ' (x)a L (x)) , 



(D2) 



(D3) 



with tf being the (adjustable) phase difference between 
the LO a,h(x) and the signal a(x). After some algebra we 
derive, on using the commutation relation (0), 
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: (X_\lu)X(lu)) = (Xj-u.')X_u< 
dxdx 1 



2tt 

[S(u) + l]I A-\ 



X{x)X{x')\e^ x '-^ 



(D4) 



where 



5 M = ^T" / dxdx ' (■X(x)X{x')) e lu{x '- x) (D5) 
2tt7 J \ I 



defines the normalized squeezing spectrum, with 



Z7T 



{al(x)d- L (x)) + (d' !i (x)a(x)) 



(D6) 



In Eq. flDq ) the symbol : : introduces normal or- 
der. We assume that the LO is prepared in a coher- 
ent state [ai,(x) — > ol(x)] and express the normally or- 
dered signal-field correlations in terms of (discrete) s- 
ordered cumulants, on recalling Eqs. (17) and ( |l8| ) and 
using Eq. ([56|). Straightforward calculation then yields 
the squeezing spectrum in the form of Eq. ( |60| ) together 
with Eqs. ( pi] ) and (|62|). 

APPENDIX E: DERIVATION OF EQ. (0) AND 
EQ. @) 

We define the operator of the number of photons in a 
sufficiently small frequency interval (ui,u! + Aui) by 



N(u) = a T (w) o(w) AloA, 



(El) 



where 



a(cu) = (27T)- 1 / 2 / dir a(x)e Ma: . (E2) 
To determine the covariance 

(AN(uj)AN(uj')) = (N(uj)N(uj')) - (N(uj)) {N(u')), (E3) 
AN(u>) = N(u>) — (N(lu)), the mean photon number 
AujA 



(N(u)) 



2tt 



dxdx' (a'(x)a(x')) 



/\\ iu(x -x) 



(E4) 



and the photon-number correlation 



{N{uj)N(d')) = 



AloA 
2vr 



x (a\x)a{x')a\y)a{y')) 



dxdx' dydy' 



!\\ iu(x' -x) iu' (y' -y) 



(E5) 



must be calculated. Using the commutation relation (^), 
Eq. (Ea) can be rewritten as 



(N(u)N(uj')) = (: N(uj)N{oj') :) + (iV(w))5(w, w'), (E6) 
where 

' 1 if \lo-uj'\ < |Aw, 



otherwise. 



(E7) 



We now express the normally ordered field correlations 
in the integrals to be performed in terms of (discrete) s- 
ordered cumulants. Recalling Eqs. ([l7|) and (|l8|), apply- 
ing the rule (pp[), and using Eq. (p6|), after some algebra 
we arrive at Eqs. ( |67| ) and (^) together with Eqs. ( |69| ) 
-(I7U). 
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